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We present a general theory of multiorbital spin waves in magnetically ordered metallic systems. 
Motivated by the itinerant magnetism of iron-based superconductors, we compare the magnetic 
excitations for two different scenarios: when the magnetic order either sets in on the on-site orbital 
level; or when it appears as an electron-hole pairing between different bands of electron and hole 
character. As an example we treat the two-orbital model for iron-based superconductors. For small 
magnetic moments the spin excitations look similar in both scenarios. Going to larger interactions 
and larger magnetic moments, the difference between both scenarios becomes striking. While in 
the excitonic scenario the spin waves form a closed structure over the entire Brillouin zone and the 
particle-hole continuum is gapped, the spin excitations in the orbital scenario can be treated as spin 
waves only in a close vicinity to the ordering momenta. The origin of this is a gapless electronic 
structure with Dirac cones which is a source of large damping. We analyze our results in connection 
with recent neutron scattering measurements and show that certain features of the orbital scenario 
with multiple order parameters can be observed experimentally. 

PACS numbers: 74.70.Xa, 75.10.Lp, 75.30.Fv, 75.25.Dk 

I. INTRODUCTION 

In recent years correlated multiorbital electronic systems have attracted a lot of attention due to a variety of 
ordering phenomena such as unconventional charge and magnetic ordering in Na^CoC^, and recently discovered 
high temperature superconductivity in iron-based superconductors^. Many of these multiorbital systems (MOS) are 
metallic and have Fermi surfaces (FS) consisting of multiple pockets, which makes an effective description with single 
band models in principle impossible. Moreover, due to strong hybridization effects, the conducting bands are usually 
linear combinations of different orbitals. It is therefore desirable to calculate spin excitations in magnetically ordered 
MOS. 

In iron-based superconductors the role of the underlying orbitals in the formation of the magnetic order and 
superconductivity is a subject of ongoing research. In particular, initial models based on the pure fermiology and the 
fact that the FS consist of hole, a, and electron, j3, pockets separated by the antiferrromagnetic (AF) wave vectors 
Qi = (7r, 0) and Q2 = (0, 7r) predicted a superconducting instability with a so-called s H —wave symmetry of the 
superconducting gap^^. Correspondingly, the spin density wave order originate s from nesting of the a and fi bands 
which yield a logarithmic instability of the spin response at the AF wave 

vectoJ^EEEo]. Within this scenario, often 
called excitonic, the role of the lattice and orbitals is negligible as the main effect is associated with the electron and 
hole character of the bands involved. Moreover, the RG studies support this indirectly as the bare local interactions 
renormalize drastically during the RG flow as a result of the interfering logarithms in the particle-hole and particle- 
particle channels^'. Note that despite its simple form, the excitonic scenario describes well spin wave excitations 
below Tn, their dispersions and anisotropy of the spin wave velocities^. 

At the same time, it was argued that the orbitals do play an important role, namely, by introducing a significant 
angular dependence for the intraband and interband interactions resulting in a change of important details in the 
superconducting and antiferromagnetic states. For example, it was argued that the originally nodeless — super- 
conducting state may acquire nodes on the electron pockets^^l. For the AF state it was predicted that iron-based 
superconductors show nodal elementary excitations with a linear Dirac spectrum at low energies as a result of orbital 
effect g^HlZI i n addition, due to the fact that the antiferromagnetic order is either (it, 0) or (0,7r), the magnetic 
anisotropy induces orbital anisotropies which in turn affect the electronic properties alon g a and b directions'" 2 ^. 
They may also introduce novel features such as ferro-orbital or antiferro-orbital ordejplBH associated with intra and 
inter-orbital electron-hole pairing. Therefore, excitonic band and orbital lattice descriptions are not equivalent even 
in the limit of interactions being much smaller than the bandwidth. Furthermore, nesting features do not play an 
essential role in the formation of the AF order starting from the orbital description. 

In this paper we analyze itinerant spin excitations of the AF state with either (it,0) or (0,7r) ordering wave vector 
for the two orbital model in the excitonic and orbital scenarios. The aim of the present work is two fold: on the one 
hand we want to present the general spin wave theory in MOS with itinerant magnetism in a self-consistent random 
phase approximation (RPA)-type formalism. We compare in detail the excitonic and orbital AF scenario and highlight 
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important differences in the properties of the magnetic excitations in both cases. On the other hand, we analyze the 
magnetic state of iron-based superconductors. In particular, we show that for small and intermediate values of the 
magnetic moment, spin waves in the excitonic scenario have longer life times while the low-energy continuum reveals 
stronger damping of them in the orbital description. We further analyze the magnetic excitations of the two orbital 
model and discuss its connection to recent experiments. In the orbital resolved susceptibilities, we observe major 
differences between the spin waves coming from different orbitals and discuss its relevance for resonant inelastic X-ray 
scattering (RIXS). 

The paper is organized in the following way: After giving the basic definitions we review the multiorbital formalism 
for calculating RPA spin susceptibilities in the paramagnetic phase and extend it to the symmetry broken AF state in 
a self-consistent way. Using as an example the two orbital modeP^we provide explicit expressions of the susceptibility 
tensor and interaction matrix and explain in detail the differences between the excitonic and orbital AF scenarios. 
Although we treat a concrete two orbital model, the formalism is easily applicable to different MOS. In the Results 
section we present a comparison between the orbital and excitonic spin waves, their damping, and evolution from 
small to large magnetic moments. Finally, we show the orbital resolved susceptibilities and discuss the possibility of 
observing them in RIXS experiments and make connection to recent INS experiments. We close the paper with a 
Summary section. 

Our results for iron-based superconductors are complementary to previous studies. For example, Ref l25l discusses 
spin waves in iron-based superconductors using the orbital picture within a five orbital model. The collective spin 
excitations in SDW phases using simplified band models have been calculated in Refs l9ll(?l Here, we analyze both 
orbital and band descriptions within the same two orbital model and critically compare both scenarios. We also 
provide an extensive discussion of the method which can easily be implemented for other multiorbital models. 

II. THEORY 

We begin by discussing general aspect of the self-consistent RPA-type theory of spin excitations in the AF ordered 
phase. For a single band model the formalism was presented in Refs.26 27] Here, we demonstrate the theory in 
detail for a two orbital model in the excitonic and orbital AF scenarios using the model Hamiltonain proposed by 
Raghu ct alPsl for iron based superconductors. The generalization to different MOS and models with more orbitals is 
st raight forward . 

A. Basic Definitions 

We define the spin operator as 

S,(q) = J2 C(k + q)%^ CT <(k), (1) 

where ip\<j is the creation operator of the orbital s with spin a in the MOS. The transverse spin susceptibility can be 
calculated from the spin-spin correlation function in Matsubara frequency space: 

rP 

xi(q,iO) = / dre^(T T 5+(q,r)5 t -(-q,0)>. (2) 
Jo 

In MOS it is convenient to introduce a dynamical spin susceptibility tensor in a more general way using four orbital 
indices 

X± gt (q,ift) = f dre** £ (T T ^ 7 (k,r)^(k+ q,r)^ y (k', 0)^(k' - q,0)) a+ a~ l5 ,. (3) 

k,k', 7,(5,7', <5' 

With the definition of the Green's function (GF) 

G° t (k,icj n ) = - J Q P dre iUT (r T Vw(k,r)v4 CT (k,0)} (4) 
the susceptibility matrix can be written as 

X^,t(q,ffi) = -j^n Gl P (K^n)G l qt {k + q,iuj n +in). (5) 
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The susceptibility in the paramagnetic phase within the RPA can be written in a matrix form see for example 
RefH3] 

RPA _ n , ~,RPA TT ..0 f R \ 

Xbsta — Xbsta'T Xb'sta' u c'b'a'd'Xbc'd'a > \°) 

and the solution of the RPA in matrix form is straightforward. 

IX} RPA = X°(1-UX )- 1 (7) 
The physical susceptibility which is measured in experiment is then obtained 

x^M = \H^ t A M- (8) 

St 

Generally, for MOS, the kinetic or tight-binding part Hq of the Hamiltonian H = Hq + Hi nt , including the intra- 
and inter-orbital hoppings, can be diagonalized via a unitary transformations which gives the eigenergies E u (k) 
of the system. The operators in the orbital and the band basis are connected via the coefficients of the unitary 
transformations a s u {k,a) 

^Mk) = ^<(k,<7) 7l -(k) (9) 
Without interactions the GFs are easily obtained as 

^) (10) 

and the transverse spin susceptibility becomes: 

± , , s 1 \- <(k, tH* (k, t)<(k + q, |R* (k + q, |) 

Lp,ilq ' 1 N ?hh [--^(k)n-+^-^(k+q)] ■ 

The Matsubara sum can be performed via a standard contour integration and by analytic continuation ifl — > f2 + iT. 
Observe that for large Hamiltonian matrices the coefficients a„ have to be calculated numerically by diagonalizing 
the Hamiltonian matrix written in k-space. In the following we move to the AF state. 

In a mean-field treatment of the SDW phase, the interaction part of the hamiltonian Hi nt is decoupled with respect 
to a groundstate with a finite magnetization at wave vector Q, e.g. (ir,0) or (0, w) for iron-based superconductors. 
The order parameters read 

aA uv , = £(VL(k)<Wk+Q)) (12) 

k 

where a = ±1 refers to up or down spins. In the magnetic state, SU(2) spin symmetry is broken which yields 
X ± 7^ \x zz ■ I n addition, the AF state also breaks translational symmetry and the mean-field Hamiltonian in the 
reduced BZ takes the form 

H SDW = £ ^t (k) \ V 3 CT ( k) (13) 



where we introduce the new spinor 



- Ut%)) ■ < i4) 

Note for MOS each component of the spinor may acquire a spinor structure in a orbital space as well. Then the GF 
acquires off-diagonal elements in the momentum space and becomes a matrix. In short-hand notation, we write only 
the orbital index with a tilde to indicate that the momentum is shifted by the AFM wave vector. We thus obtain 
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Analogously, the susceptibility matrix becomes a matrix in q-spao 



Xspqt 



( xf pq t (q» q) xfp q t (q, q + Q) 
vxf pgt (q + Q,q) xt Mt (q + Q>q + Q) 

G S pGqt + G%%G 



GspGqt 

GspG q i ■ 



GspGqt 



T sp^qt 
GspGqt 



GspG q t 
GspGqt 



GspGqt 



GspGqt 



GspGqi 4- G S pG qt 



GspGqt + GspGqt 
GspGqt + GspGqt 



(16) 
(17) 



where the off-diagonal components are called Umklapp susceptibilities as they connect spins with momentum shifted 
by the AFM wave vector, which is a reciprocal wave vector of the magnetic BZ. 

The form of the RPA in the SDW state remains the same, see Eq.([6]), with the only extension that the interaction 
matrix also doubles in size with U matrices on the diagonal and zeros on the off-diagonal. To obtain the physical 
susceptibility we perform an RPA with Eq.(|17j) and trace over the orbital indices of the first entry xl^qf (^i q) 01 
the susceptibility matrix. The imaginary part is directly related to the cross section of INS measurements. We also 
recall that self-consistency enters via taking the same set of interaction parameters in the RPA and the Hartree-Fock 
mean-field equations. 



B. Two orbital model: excitonic versus orbital SDW scenario 



Let us illustrate the derivation, presented above, for the two orbital modeP^with d xz and d yz orbitals proposed for 
the iron-based superconductors. The kinetic part of the Hamiltonian is written in the two component spinor basis in 
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orbital space with the help of the Pauli matrices r. 

Ho = ( k ) K £ + ( k ) - to + e_ (k)r 3 + 6,, (k) n ] (k) 



k.rr 



where 



Vv(k) 



dxza (k) 

dyza (k) 



(18) 



(19) 



Here, e±(k) — g^WjEgvW ; £ x (fc) = — 2fiCosfc. r — 2t 2 cos/c y — 4i 3 cos fc^ cos fc y , e y (fc) = — 2i 2 cosfc. r — 2ticosfc J/ 
4i 3 cos fca; cos fc y , and £ xy (k) = — 4i4 sin fc^ sin The resulting Fermi surface is shown in the inset of Figflj 
The interacting part of the Hamiltonian consists of different on-site electron-electron interactions: 



H Zt = ' ' 1 X! n "" T " 

1 V V^fJ,,<T,<T' 



11/(7 ' HfJUf" 



(20) 



where C/ and V refer to the intra- and inter-orbital Coulomb repulsion, J and J' = J/2 denote the Hund and the pair 
hopping terms. The susceptibility matrix, with labels 1 for xz and 2 for yz, reads 



Xini X2111 X1112 X2112 

X1211 X2211 X1212 X2212 
X1121 X2121 X1122 X2122 
.X1221 X2221 X1222 X2222. 

and the interaction matrix for the transverse component of the spin susceptibility is given by: 



U = 



(21) 



-u 








J- 





J 


U-2J 








U-2J 


J 





J 








17. 



(22) 



Following the procedure described above it is straightforward to compute the physical susceptibility in the paramag- 
netic state. In particular, Figjl] shows the real part of the physical susceptibility for the two orbital model. Observe 
that in the paramagnetic state the response is governed by the antiferromagnetic fluctuations centered at (it, 0) and 

(0,7T). 

In the next step we introduce self-consistent mean-field averages associated with the AF state. As described above, 
we consider two scenarios. 

In the first case, which we denote as orbital, one introduces the magnetic state in the orbi tal ba sis using the 
local interactions. The mean-field phase diagram of the two orbital model was studied extensivel) * 15 * 21 ! and the main 
conclusion is that AF ordering involves only very weakly interorbital components. Thus we assume in the following 
that the mean-field AF order involves only states within the same orbital, separated by Qa_f (&vv> — <^/i/'A„), and 
decouple the on-site local interaction, Eq.(20), with respect to it. In addition, we set Qaf = (0, 7r) and do not consider 
the problem of its selection out of the manifold of degenerate magnetic states as it was considered by us previously^. 
We obtain two coupled gap equations 



M x = -[UA XZ + JA yz ] 
M y = — [UA yz + J A xz ] 



and the Hamiltonian becomes 



' Orbital 







crM x 







Ey 





aMy 


<jM x 


Q 









OMy 


£xy 


Ey 



(23) 
(24) 
(25) 



(26) 



where the tilde again means that the momentum is shifted by Qaf The matrix can be inverted and with the new 
eigenvectors one can solve the coupled gap equations and carry out the spin wave calculations as outlined above. 

The second, excitonic, scenario is built on the belief that in the iron-based superconductors, the antiferromagnetic 
instability, at least partially, is connected to nesting of electron and hole bands^l. In this scenario, the log instability 
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FIG. 2: Fermi surface of the AF state with a magnetic moment of m = 0.23/is in the orbital (a) and excitonic (b) scenarios. 
Despite the fact that the FS topology is similar in both cases its origin is different. While in the excitonic scenario the remnant 
FS pockets arise due to imperfect nesting between electron and hole bands and a small size of the SDW gap, the electronic 
structure and the Fermi surface pockets in the orbital scenario are consequences of the nodal and the pockets do not 

vanish for increasing the size of the SDW gaps but simply shift their position in the reduced BZ. 



of the particle-hole response function is similar to that in the particle-particle Cooper-channel. In this regard it 
seems natural to assume that the magnetic instability arises not in the orbital representation but on the band level for 
electrons lying close to the Fermi level. This situation is not unique to ferropnictides and is believed to arise also in Cr 
and is sometimes dubbed excitonic antiferromagnetism. Within this scenario the tight-binding hamiltonian,Eq.( 18 1, 
is first diagonalized to obtain the band structure. One further assumes, that the bands are the relevant degrees of 
freedom and ignores the role of the underlying orbital structure. The interactions are then also re-written in the band 
representation. Moreover, within the RG flow they are strongly renormalized, to such an extent that their values at 
the Fermi energy are not the same as the initial ones determined by the local on-site interactions^^. This allows us 
to approximate the band interactions as Fermi-liquid density-density type interactions, which at the lowest order are 
angular independent. Following this we consider the density-density interactions which give rise to SDW ordering at 
(7T,0) or (0,tt^H 



H 



band 



E 

k,k',q 

cr, cr' 



"1 



7i, 



_(k + q) 7 ^(k' 



q)72<r'(k')7i<r(k) 



+ i 



7 L(k + q) 7 L(k' - q)72a'(k')72 ff (k) + h. c .(27) 



The 7-operators, see Eq.([9]), are written in the band representation which diagonalizes the tight-binding part of the 
Hamiltonian with Eigenenergies -E„(k). 

In the excitonic scenario the AF order involves quasiparticle states from the electron and hole bands connected 
by the AF momentum Qx = (0,7r). By decoupling the band interaction with respect to A 12 one finds a single gap 
equation 



M = - 



u x + u 3 



(28) 
(29) 



and the Hamiltonian matrix becomes 



/ E 1 aM\ 

E 2 oM 

oM E x 

\crM E 2 J 



(30) 



Within the excitonic scenario we obtain a constant gap splitting of the excitation spectra all over the first BZ. By 
contrast, in the orbital picture the overall gap in the excitation spectra appears to be nodal as a result of the underlying 
symmetry of the orbital structure^. In addition, the structure of the AF state is richer as it involves at least two 
intra-orbital SDW orders. Despite these differences between the orbital and excitonic scenario, it is interesting to 
note that the resulting FSs for small sizes of the magnetic moments appear to be quite similar as shown in Fig(2] 
Nevertheless, the origin of the pockets is different. In the excitonic scenario the remnant FS pockets arise due to 
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FIG. 3: Calculated sizes of the magnetic moments in the excitonic and orbital scenarios as a function of the interaction strength. 
The solid and the dashed lines are guides to the eye. 



imperfect nesting between electron and hole bands and a small size of the SDW gap. On the contrary, the Fermi 
surface pockets in the orbital scenario are a consequence of the nodal SDW^ and they do not vanish for increasing size 
of the SDW gaps but simply shift their position in the reduced BZ. The difference between both scenarios becomes 
apparent with increasing size of the magnetic moments as, we show below. 

In Fig|3]we show the evolution of the magnetic moment within the orbital and excitonic scenario with increasing 
interaction strength (U + J and u\ + u% respectively). Observe that the excitonic scenario requires a smaller value of 
the interaction to obtain sizeable SDW gaps and the corresponding magnetic moments. This is a result of neglecting 
orbital matrix elements which overall tend to weaken the effect of interactions in the SDW gap equation. Indeed within 
an excitonic scenario, the SDW gap and the magnetic moment appear for smaller magnitude of the interactions as 
compared to the orbital case. The origin is the logarithm of the particle-hole channel in the excitonic scenario. 
Although, also present in the orbital scenario, its effect is weakened by the underlying orbital character of the wave 
functions. Thus within the orbital scenario the SDW gaps increase much slower with increase of the moment. In 
addition for /i > 0.38/is the Fermi surface is completely gapped within the excitonic scenario, while FS pockets are 
always present in the orbital scenario. The latter leads to finite damping of the spin waves; well-defined dispersive 
spin waves are thus absent from much of the BZ in the orbital scenario as will be shown in the following. 



III. RESULTS 



In the following we analyze spin excitations within the excitonic and orbital scenarios. The interaction parameters 
- U and J for the orbital scenario, and u\ + for the excitonic scenario - determine the magnitude of the magnetic 
moment via the self-consistent gap equations, Eqs.(28 and 23). We always choose them in such a way that the 
resulting moment is equal within both scenarios which allows a direct comparison of the spin waves. We concentrate 
on the undoped case that corresponds to a total number of particles n — 2 per unit cell. Observe also that in the 
orbital scenario occupation and magnetization of the two orbitals are unequal because the AFM order with Qi=(0, tt) 
ordering momentum breaks the symmetry between the two orbitals. 

In Fig|4]we show the evolution of the transverse spin excitations in the orbital (left) and excitonic (right) scenarios 
along the route (tt, 0) -+ (0,0) -+ (0, tt) — > (tt, tt) of the first BZ and for various values of the magnetic moment. The 
intensity is shown on a log intensity scale. In panels (a)-(b) the magnetic moment is small, m = 0.23/is. In the 
excitonic case we use \i = 1.54i, (u\ + u%)/2 = 1.335i which gives a gap of A12 = 0.1525i, while in the orbital case we 
employ U — 3.5i and the gaps are A xz = 0.2086i and A yz = 0.2283i. As mentioned above, larger magnitudes of the 
gaps are required in the orbital case. 

The peculiar feature of the spin excitations in the AF state within both scenarios are the well-defined spin wave 
excitations around the ordering wave vector Qi = (0,7r) and their visible dispersion up to energies of about O.li 
despite the finite damping originating from remnant FS pockets. At higher energies the spin waves are overdamped by 
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FIG. 4: Calculation of the imaginary part of the transverse spin susceptibility in the orbital (left panel) and excitonic (right 
panel) scenarios for the magnetic moment of m = 0.23/is(a)-(b), m — 0.38^b(c)-(c!) and m = 0.75/is(e)-(f). The intensity in 
states/eV is shown on a log scale. For concreteness, we set J = 0.1(7 and use a quasiparticle damping of T =0.005t. 

the particle-hole continuum and evolve further as Stoner-like excitations (paramagnons). Another common element 
of both scenarios for this small value of the magnetic moment is that the particle-hole continuum dominates the 
low-energy excitations in the entire BZ except around the ordering wave vector itself. Although the details of the 
continuum slightly differ in both scenarios, the gross features are the same. In particular, the low-energy excitations 
around (0, tt/2), (tt, 7t), and (n, 0) are due to particle-hole scattering between remnant parts of the Fermi surface, see 
Fig(2] The higher energy features look slightly different due to the richer structure of the particle-hole continuum in 
the orbital scenario, but their intensity is overall low. 

The difference between both scenarios becomes more apparent if one raises the moment towards m = 0.38/is, see 
FigHtc)-(d). Here, we use fj, = 1.53t and (itj + u 3 )/2 = 1.465i with the corresponding SDW gap of A 12 = 0.2803i, 
while in the other case one needs U = 3.95t and A xz = 0.4008£ and A yz = 0.4250t. Observe that in this case the 
continuum of the spin excitations is completely gapped in the excitonic scenario, as shown also in Fig[5jb). This is 
a result of the fully gapped Fermi surface in this case. As a consequence, the spin wave branch is well-defined over 
almost the entire range of the calculated momenta. 

It is interesting to note that within the two orbital model and excitonic scenario of the SDW there is an 0(6) 
degeneracy of the magnetic ground state^ which manifests itself in additional zero energy modes at the non-ordered 
momenta of (tt, 0) and (tt, tt) as shown in FigQd). These modes, however, are accidental in the sense that they are 
not stable with regard to perturbations such as the finite Landau damping introduced by the remnant Fermi surfaces. 
This is why they are not seen for m — 0.28^#. However, as soon as the FSs are completely gapped they become 
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visible. In the orbital scenario, however, the FSs will never be gapped completely for symmetry reasons and the 
particle-hole continuum remains gapless as shown in Figj5^a). Thus, the spin waves remain well-defined only around 
the ordering momentum. Away from it they are overdamped because the Dirac cones in the electronic excitation 
spectra lead to remnant FS pockets. 

The situation does not change significantly upon further increase of the interaction. In particular, in Fig[4|e)-(f) 
we show the results for the transverse part of the spin susceptibility for the magnetic moment of to = 0.75/is. This 
magnetic moment arises for /i = 1.53£ and (u\ + u 3 )/2 = 2.043/j with A 12 = 0.7607i in the excitonic scenario and 
for U — 5.07;, A xz = 0.9846£ and A yz = 1.0772/; in the orbital ones. The essential features of the spin waves in the 
orbital and exciton scenario remain the same. While in the orbital picture the spin waves are overdamped away from 
the ordering momentum, the spin waves in the excitonic scenario form a closed structure over the entire BZ. Observe 
also that within both scenarios the spin waves show an anisotropy along the antiferromagnetic (v y ) and ferromagnetic 
(v x ) directions. Namely we we find that within both scenarios the anisotropy is relatively weak, v y Jv x « 1.04 in the 
excitonic and v y /v x « 1.06 in the orbital scenario. This agrees with the result of the J\ — J2-model2ZI. 




(n,0) (0,0) (0,71) (ivi) "(™»0) (0,0) (0,n) (7t,7t) 



FIG. 5: Imaginary part of the bare transverse physical susceptibility, Imx°, for m = 0.38/iB for the orbital (a) and excitonic 
(b) scenario with m = 0.38/i_g- 

A peculiarity of the spin waves within the orbital scenario is the relatively rich structure of the spin excitations 
away from the ordering momentum. Here the complication is that each orbital has its own magnetic gap which 
correspondingly affects the magnetic excitations within each orbital. In Figj6] we display the intra and inter-orbital 
components of the transverse spin susceptibility for the orbital and the excitonic scenario of the SDW phase on 
an absolute intensity scale. Observe that within the orbital scenario, the Goldstone modes originate from the intra- 
orbital xz and yz susceptibilities, while in the excitonic scenario the entire low energy excitations are in the inter-band 
susceptibility by construction. In addition, as one sees from the left panel, (0, it) order breaks the symmetry between 
xz and yz orbitals which leads to different behavior of the spin excitations away from the ordering momentum and 
higher energies visible in the corresponding intra-orbital susceptibilities. Given the fact that the spin waves show 
some orbital dependences we believe that this feature could be measur ed by RIXS experiments^ 21 . 

Another interesting feature of the spin excitations in the orbital scenario is that their spin wave branches reappear 
out of the continuum at higher energies, e.g. when going from (0,0) to (0,tt) as shown in Figj4^c),(e). Here, the 
orbital coherence factors yielding the larger continuum (damping) are small at these higher energies, thus allowing 
the spin waves to reappear. As soon as the itinerant magnetism appears in MOS with multiple order parameters this 
reappearance of the spin waves at higher energies away from the ordering momentum should be a general feature not 
restricted only to iron-based superconductors. 

Spin excitations within the orbital scenario show a rich structure at energies which exceed well the onset of the 
particle-hole continuum. Note that at these energies the spin excitations should not any longer be interpreted as spin 
waves originating from breaking of the spin-rotational symmetry in the antiferromagnetic state but rather should be 
attributed to a Stoner continuum. These excitations are almost unchanged across the Neel temperature and as such 
can be well separated in the experimental data. 

In Figj7]wc show the behavior of the Im^ summed over the q x and q y momenta (compare to FigJH]) for two different 
magnetic moments. Observe that for to = 0.22/jb (left) the damping increases slowly at low frequency and remains 
largely featureless up to energies of u> ~ 0.08t. At higher frequencies damping becomes significantly larger and all 
excitations above this frequency refer to the Stoner continuum. For m — 0.38/is (right) the continuum is gapped in 
the excitonic scenario at low frequencies but apart from this, the situation remains the same. In particular, at low 
enough frequencies, damping is weak and spin excitations can be associated with spin waves. For energies higher 
than u ~ 0.4?; the damping is so large that the excitations exist in the form of the Stoner continuum only. The 
latter will only be weakly sensitive to the antiferromagnetic transition. Note that several experiments have studied 
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FIG. 6: Components of the transverse RPA spin susceptibility for the orbital (left panel) and excitonic (right panel) scenarios 
for a magnetic moment of m — 0.38/j.b- The intensity is shown on an absolute scale. 
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the temperature dependenc e! 38 * 39 * of the spin excitations and found that the high energy part of the spectrum is not 
sensitive to the magnetic or structural transition. In contrast to their interpretation of this result as local moments 
present at higher energies, we believe that the most natural explanation of these experiments should be given in terms 
of the Stoner continuum which despite its structure does not refer to spin wave excitations originating from the AF 
state but is mostly a property of the paramagnetic state. Therefore, it is almost unchanged across the transition. 

Despite the widely known fact that the two orbital model is too simplistic to lead to quantitative results there arc 
several aspects of our calculations which we believe will also exist in a more accurate five orbital model. This includes 
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the anisotropies of the spin wave velocitie s] 34 * 35 !, occurrence of strong damping of the spin excitations in the orbital 
scenario with multiple order parameters and the hidden symmetries of the spin waves within the excitonic scenario 
which manifest themselves in accidental zero modes at certain non-ordering momenta like (tt,tt) or (tt, 0). We also 
find that for the nodal SDW, it is unavoidable to have a particle-hole continuum even at low energies away from the 
ordering wave vector independent of the particular details of the model. The puzzle wether or not damping of spin 
waves via the particle-hole continuum is observed can be resolved by the observation that in our calculations damping 
is always there even at low energies as long as some remnant portions of the FS are present. As seen in FigjTJ the 
continuum itself can have additional steps and structure as a function of frequency. 

Regarding whether the pure excitonic or orbital scenario has to be employed depends on the energies. For example, 
recently Ewings et alP^ measured spin waves in SrFe2As2 throughout the whole BZ and up to high eneries. They 
conculded that a purely local moment model is not able to fit the data. By comparing excitonic scenario^ with 
the orbital one^ they found that the low energy part of the spin wave dispersion can be well described by both 
scenarios, whereas, the orbital scenario gives generally a better description of the spin excitations at high energies. 
Most likely an inclusion of the momentum dependence of the interactions in the excitonic scenario would also improve 
the comparison. 



IV. SUMMARY 



Motivated by the magnetism present in the metallic iron-based superconductors we have analyzed spin excitations 
in magnetically ordered multiorbital systems within a self-consistent mean-field RPA-type formalism. We presented 
an extensive discussion of this theory and as an example carried calculations for the two orbital model. This formalism 
can be easily generalized to different MOS. As the magnetic electron- hole pairing can either reside on the orbital or 
at the band level, we have compared the behavior of spin excitations in two scenarios: either the on-site orbital one, 
or the excitonic scenario, where the magnetic order occurs on the band level in momentum space. For small magnetic 
moments the low energy spin excitations look similar in both cases. In particular, the excitations are spin waves 
around the ordering momentum which become overdamped at higher energies away from Qaf- Due to remnant 
FS pockets the behavior of the particle-hole continuum is quite similar in both cases. The situation changes for 
larger magnetic moments, as the particle-hole continuum gets completely gapped within the excitonic scenario, while 
it remains by symmetry reasons gapless due to Dirac cones in the orbital description. As a result in the excitonic 
scenario spin waves develop over the entire BZ with accidental zero modes at (jr, 0) and (jr, it) points of the BZ due 
to an underlying accidental 0(6) symmetry of the spin excitations in this case. At the same time, in the orbital 
picture the spin waves remain well-defined only around the ordering momentum and are overdamped away from it. 
In addition we also find that due to multiple gaps in the orbital scenario the spin wave branches may reappear deeply 
inside the continuum because of the orbital character of the wave functions. Analyzing the components of the spin 
susceptibility and the behavior of the particlc-hole continuum within both scenarios one could clearly distinguish the 
orbital and the excitonic origin of the spin excitations in iron-based superconductors. 

We thank A.V. Chubukov, P.M.R. Brydon, C. Timm, P. Dai, and T. Perring for useful discussions. The research 
of JK was supported by a Ph.D. scholarship from the Studienstiftung des deutschen Volkes. IE acknowledges the 
DAAD PPP Grant N50750339 for the financial support and MPI PKS for hospitality. 
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